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ABSTRACT 

We have reanalyzed data from observations of PSR B1706-44, SN 1006, and the Vela pulsar region 
made with the CANGAROO 3.8 m Imaging Atmospheric Cherenkov Telescope between 1993 and 
1998 in response to the results reported for these sources by the H.E.S.S. Collaboration. Although 
detections of TeV gamma-ray emission from these sources were claimed by CANGAROO more than 
10 years ago, upper limits to the TeV gamma-ray signals from PSR B1706-44 and SN 1006 derived 
by H.E.S.S. are about an order of magnitude lower. The H.E.S.S. group detected strong diffuse TeV 
gamma-ray emission from Vela but with a morphology differing from the CANGAROO result. In our 
reanalysis, in which gamma-ray selection criteria have been determined exclusively using gamma-ray 
simulations and OFF-source data as background samples, no significant TeV gamma-ray signals have 
been detected from compact regions around PSR B1706-44 or within the northeast rim of SN 1006. 
The upper limits to the integral gamma-ray fluxes at the 95% confidence level have been estimated 
for the 1993 data of PSR B1706-44 to be F(> 3.2 ± 1.6TeV) < 8.03 x lO-^^ photons cm-^s"!, 
for the 1996 and 1997 data of SN 1006 to be F{> 3.0 ± 1.5 TeV) < 1.20 x lO^^^ photons cm^^ g-i 
and F{> 1.8 ± 0.9 TeV) < 1.96 x 10~^^ photons cm~^s~^, respectively. We discuss reasons why 
the original analyses gave the source detections. The reanalysis did result in a TeV gamma-ray 
signal from the Vela pulsar region at the 4.5cr level using 1993, 1994, and 1995 data. The excess 
was located at the same position, 0?13 to the southeast of the Vela pulsar, as that reported in 
the original analysis. We have investigated the effect of the acceptance distribution in the field of 
view of the 3.8 m telescope, which rapidly decreases toward the edge of the field of the camera, on 
the detected gamma-ray morphology. The expected excess distribution for the 3.8 m telescope has 
been obtained by reweighting the distribution of HESS J0835-455 measured by H.E.S.S. with the 
acceptance of the 3.8 m telescope. The result is morphologically comparable to the CANGAROO 
excess distribution, although the profile of the acceptance-reweighted H.E.S.S. distribution is more 
diffuse than that of CANGAROO. The integral gamma-ray flux from HESS J0835-455 has been 
estimated for the same region as deflned by H.E.S.S. from the 1993-1995 data of CANGAROO to be 
F{> 4.0 ± 1.6 TeV) = (3.28 ± 0.92) x lO^^^ photons cm^^g-i^ .^^j^ich is statistically consistent with 
the integral flux obtained by H.E.S.S. 

Subject headings: gamma rays: observations — methods: data analysis — pulsars: individual 
(PSR B1706-44, Vela pulsar) — supernova remnants — supernovae: individual 
(SN 1006) 
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The CANGAROO 12 3.8 m telescope was an Imag- 
ing Atmospheric Cherenkov Telescope (lACT) for TeV 
gamma-ray astrophysics which operated from 1992 to 
1998 near Woomera, South Australia. It is retro- 
spectively called the CANGAROO-I telescope since 
its pioneering role as the first lACT in the south- 
ern hemisphere was subsequently inherited by the 
CANGAROO-II telescope (TanimOTi et al. 19 9j) and th e 
CANGAROO-III telescope system (IKubo et al1[200l . 
The discoveries of TeV gamma-ray signals from five 
sources have be en claimed from C ANGAROO-I data: 
PSR B 1706-44 (iKifune et a.l.l fl99l. the Vela puls ar re- 
gion (|Yoshikoshi et al.l ll997[), SN 1006 jT animori 'et all 
1998), RX J1713J:J946_|Muraishi_it3j [2000), and 
PSR B1509-58" (jSako et aL 2000) . However, H.E.S.S. 
has later reported upper limits to the TeV gamma-ray 

Acronym for the "Collaboration of Australia and Nippon 
(Japan) for a GAmma Ray Observatory in the Outback" . 

A marginal detection compared to the other CANGAROO-I 
sources. 
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fluxes fro m PSR B 1 706-44 (lAharonian et al.ll2005aD and 
SN 1006 (jAharoiiian et al.ll2005H ). which are lower than 
the CANGAROO-I fluxes by factors of ~ 10. The re- 
sults are obviously inconsistent as the reported statisti- 
cal significances of the CANGAROO-I detections were 
12cr for PSR B1706-44 and T.Tcr for SN 1006, assum- 
ing that the flux levels are constant over the 10 year pe- 
riod. More recently, H.E.S.S. has detected a TeV gamma- 
ray signal from SN 1006 aft er accumulating more data 
(|Naumann-God6 et al.l 120081. but the flux level is well 
below their previous upper limit. The Durham group 
also reported a statistically significant (5 .9g) detection of 
TeV gamma rays from PSR B1706-44 (jChadwick et al.l 
|1998() . Preliminary support of the TeV gamma-ray sig- 
nals from PSR B1706-44 and SN 1006 were reported 
from CANGAROO-II observations (K ushida et al.ll2001l . 
[200a iHara eFaI][200lh . 

H.E.S.S. also observed the Vela pulsar region in 2004 
and 2005, and detected a stron g gamma-ray signal at th e 
50% Crab level above ITeV (|Aharonian et all l2006al ). 
The emission is diffuse over the X-ray (0.9-2.0 keV) 
j et-like image, which w as first detected by ROSAT 
(|Markwardt &: OgelmanI 119951 ). and extends about 1° 
from the pulsar toward the south-southwest, reaching 
the position of the brightest radio emission i n the Vela 
super nova remnant (SNR), Vela X ( Wei ler fc Panagial 
Il980l) . The position of the maximum TeV emission 
in HESS J0835-455 is located in the middle of the 
jet-like feature, which has been identified by Chandra 
to be ejected in the direction of the equator ial pulsar 
wind (,Hclfan d et all 120011 : iPavlov et al] 120031 ). On the 
other hand, the TeV emission detected by CANGAROO- 
I wa s offset from the pul sar by 0?13 to the south- 
east (jYoshikoshi et al.lll997D . corresponding to the pos- 
sible pulsar birthplace calculated using the proper mo- 
t ion of 45 mas yr~ ^ toward a position angle of 301° 
Dodson et al.ll2003|) and the characteristic age of 11 kyr 
Manchester et al.l l2005f) . The CANGAROO-I source 



thus appears to be different from the H.E.S.S. source, as 
the angular offset between their peaks is 0?34, although 
the position of the CANGAROO-I peak does lie within 
the region of the H.E.S.S. diffuse emission. 

The CANGAROO-I resufrs for Vela obtained by 
I Yoshikoshi et al.l (|1997f ) were based on the data taken 
in 1993, 1994, and 1995. The Vela pulsar region was 
also observed in 1997 with the recoated 3.8 m reflec- 
tor and a TeV ga mma-ray sig nal was again detected, 
at the 4.1(7 level (lYoshikoshil I1998D. Th e 1997 data 
were later also analvzed bv iDazelev et al.l (|2001l), who 
reproduced the same sig nal as 'Yoshikoshil ()1998 l) but 
cast doubt upon its gamma-ray-likeness. No signal 
was found within a 2° field of view centered on the 
pulsar position using gamma-ray selection criteria op- 
timized a priori with the ir Monte Carlo simulations 
(jPazelev fc Patterson|[200l[) . 

In response to these inconsistencies, we have reana- 
lyzed CANGAROO-I data of PSR B1706-44, SN 1006, 
and Vela, including the data used in the previous analy- 
ses, with our best current knowledge. The following facts 
were not known or not fully considered in the previous 
analyses: (1) the effects of some noise produced by the 
electronic system were only discovered after the previous 

^•^ |http://«w«.mpi-hd.mpg.de/hfm/HESS/pages/home/som/2008/687| 



publications; (2) details of the analyses of CANGAROO- 
I data had varied from source to source, and sometimes 
from data set to data set, without a common, consistent 
treatment based on the experience gained over the life- 
time of the telescope; (3) improved calibration methods 
had not been applied to some old data; (4) although the 
fact that the acceptance for gamma rays is not uniform 
across the field of view was well understood, the effect 
on gamma-ray morphologies had not been investigated in 
detail; and (5) for Vela, a more precise energy spectrum, 
measu red by H.E.S.S., is available now (jAharonian et al.l 
l2006a^ and the systematic uncertainty of the integral flux 
has been reduced. Our aim in this paper is to investi- 
gate these issues on the CANGAROO-I data in detail, 
although the information available on the previous anal- 
yses is, in places, limited. We briefly summarize the 
CANGAROO 3.8 m lACT and its specifications related 
to the analyses in Section [2l and then the data of the 
three objects are summarized in Section[3l The common 
and standard analysis method redefined considering the 
above problems is described in Section 0] in detail, and 
the results obtained with it are summarized in Section [5l 
Several issues on the old results such as reproducibility 
and reliability are discussed in Section [SI and we finally 
conclude in Section [T] 

2. CANGAROO 3.8 M TELESCOPE 

The CANGAROO 3.8 m lACT ()Hara et al.lll993f ) was 

operated near Woomera, South Australia (136°47'E, 
31°06'S, 160m a.s.l.). The telescope detected Cherenkov 
photons from extensive air showers (EASs) generated by 
primary gamma rays and cosmic rays. The parabolic re- 
flector of the telescope had both a diameter and a focal 
length of 3.8 m. The telescope was originally used for 
lunar ranging, and the reflectivity of its Kanigen-plated 
surface was about 45% at wavelengths of atmospheric 
Cherenkov light. The reflector was recoated with alu- 
minum in 1996 October and its reflectivity improved to 
about 75%, which deteriorated to about 55% by 1998. 
The refle ctivity was monitor ed using a portable reflec- 
tometer (jPowden et al.lfr997| ) with a systematic error of 
about 5%. 

The telescope had an imaging camera of Hamamatsu 
R2248 square photomultiplier tubes (PMTs) in the fo- 
cal plane. Observations started with 224 PMTs, which 
were increased to 256 PMTs in 1995 April by adding 
34 PMTs at the corners of the camera. The number 
of PMTs decreased to 240 because of hardware trouble 
in 1998 February. Each PMT viewed an angular ex- 
tent of 0?12 X 0?12, and the field of view of the cam- 
era was about 3° across with a 0?18 spacing between 
pixel centers. Light guides were used from 1995 Novem- 
ber to capture Cherenkov photons otherwise incident on 
the dead space between photocathodes but abandoned in 
1996 November owing to complications described later. 

Outputs of the PMTs were fed into the electronics, con- 
sisting of trigger and data acquisition circuits after am- 
plification by buffer amplifiers. The data acquisition cir- 
cuits are based on analog-to-digital converters (ADCs), 
time-to-digital converters (TDCs), and scalers for indi- 

This value is the average of the whole surface and smaller than 
that of the recoated surface (~ 80%) since only the outer part of 
the reflector was recoated. 
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TABLE 1 

Summary of the CANGAROO-I data sets used in this 

PAPER. 



Target 



Year 



(hr) 



(hr) 



PSR B 1706-44 


1993=^ 


66.9 


(44.3) 


50.3 


(44.3) 




1994 


19.4 


(19.1) 


19.6 


(19.1) 




1995 


33.9 


(26.0) 


28.4 


(26.0) 




1997 


24.1 


(18.4) 


19.8 


(18.4) 




1998 


44.2 


(31.5) 


35.0 


(31.5) 


Vela 


1993'' 


52.0 


(32.6) 


42.2 


(32.6) 




1994'' 


63.8 


(41.9) 


56.7 


(41.9) 




1995*= 


50.2 


(44.9) 


46.7 


(44.9) 




1997 


28.8 


(16.6) 


21.4 


(16.6) 


SN 1006 


1996^^ 


26.8 


{ 8.2) 


14.7 


{ 8.2) 




1997=^ 


34.7 


(28.4) 


29.6 


(28.4) 



" Use d in IKifune et aU IflMB) . ^ Us ed in lYoshikoshi etaP 
HTWTh ." Used in lTanimori et al.l ||T998|) .'' Observation time for 
ON-source (OFF-source) after rejecting periods affected by clouds 
or significant electronic noise. The number in brackets is the ob- 
servation time after matching ON- and OFF-source observations 
in horizontal coordinates. 

vidual PMTs, which record brightness and arrival tim- 
ing of Cherenkov hght, and brightness of night-sky back- 
ground (NSB) hght, respectively. 

3. DATA 

The CANGAROO-I data used in this paper are listed 
in Table[TJ The data have been taken in on/off mode, in 
which OFF-source data to estimate the background level 
are taken on the same night, but with an offset in right as- 
cension, as the ON-source observation which included the 
target object in the field of view. The total observation 
times for ON- and OFF-source data in Table [T] (Ton and 
Tqff, respectively) are those after rejecting periods af- 
fected by clouds or significant electronic noise, which are 
determined from descriptions in logbooks and/or back- 
ground event rates. The ON- and OFF-source data have 
further been matched in horizontal coordinates of the ob- 
servations to make observation conditions as similar as 
possible, since some data were not taken following the 
above policy, and the corresponding observation times 
are given in brackets. In the following analyses, we use 
the data sets after matching ON- and OFF-source expo- 
sures unless otherwise specified. 

The data used in the previous papers are indicated 
by the superscript marks in the year column in Table [TJ 
PSR B1706-44 was also obse rved in 1992, and tho se data 
were used in the analysis of iKifune et al.l (|1995[ ). How- 
ever, they are not used here because some of the original 
data files have been lost and some necessary calibration 
data were not taken at that time, and it is expected that 
the performance over this first observing season will be 
more variable as problems were identified and addressed. 

4. ANALYSES 

Except where otherwise specified, a common analysis 
method is used for all data sets in this paper, since the 
nature of Cherenkov light images of EASs is independent 
of the observed object. In the previous CANGAROO- 
I papers, the analysis conditions differed for each ob- 
ject without a single, consistent underlying philosophy, 
in part as the understanding of the telescope improved 
with time. It is possible, however, to enhance an appar- 
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Fig. 1. — Calibrated tracking centers in the camera in the case 
of the Vela observations over five years. The origin is the center of 
the camera, and the Y-axis points to the zenith. Tracking centers 
of the 1993, 1994, 1995, and 1997 data are represented by solid 
circles, open circles, solid squares, and open squares, respectively. 

ent signal by using some of the many degrees of freedom 
in reducing the background level. To reduce the degrees 
of freedom in the analyses and ensure the results are ro- 
bust, we use here a simple analysis method similar to the 
"Sup ercuts" method of the Whipple group (| Punch et al.l 
Il991h . which is described in Section [44l 

4.1. Calibration of Tracking Centers 

The telescope tracking position in the field of the imag- 
ing camera is calibrated using the time variation of scaler 
counts of individual PMTs. When a bright star image is 
within the field of view of a PMT, its count rate increases. 
Since the telescope has an alt-azimuth mount, bright star 
images revolve around the center of the field of view of 
the camera, and the time profile of scaler counts shows 
peaks due to the passages of stars. The offset angle of the 
tracking direction from the center of the camera can be 
estimated by comparing the observed time profiles with 
simulated ones. It is necessary for this calibration to be 
done night by night as issues with the mount were found 
to sometimes result in offsets of greater than 0?1, com- 
parable to the size of the point-spread function (PSF). 
Figure [U shows an example of tracking centers calibrated 
night by night in the case of the Vela observations over 
five years, in which the maximum difference between the 
vertical offsets is more tha n 0?4. Deta i ls of t he calibra- 
tion procedure are given in lYoshikoshil (fl996l ). 

4.2. Image Generation and Cleaning 

The relative Cherenkov light signal of the ith pixel Si 
is obtained using the following formula: 



9i 



(1) 



where ai,pi, and g.i are the ADC value, the ADC pedestal 
value, and the relative gain of the ith pixel, respectively. 
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The ADC pedestal value of each pixel is determined us- 
ing Cherenkov light triggered events. After omitting any 
pixel which contains a signal or is adjacent to a pixel 
containing a signal, the ADC distribution of each pixel 
is fitted by a Gaussian distribution. The presence of a 
TDC signal is used to determine that a pixel contains a 
signal. The ADC pedestal value, taken to be the mean 
value of the Gaussian distribution, is subtracted from 
the ADC value. The pedestal standard deviation repre- 
sents the noise level of the ADC value due to NSB light 
and electronic noise. Any PMT signal lower than the 
IfT threshold based on the noise level is rejected from 
the image. Pixel gains relative to their average are de- 
termined using LED data, which were taken by artifi- 
cial triggers with a LED, permanently positioned at the 
front of the camera for flashing uniform light at the all 
PMTs, before or after Cherenkov light observation runs. 
Cherenkov light images are flat-fielded, i.e., divided by 
Qi as in equation ([T]). 

Background levels of hadronic events estimated us- 
ing OFF-source data are biased by the difference of 
NSB brightnes s between ON - and OFF-sources. Soft- 
ware padding (|Cawlevlll993f l is an effective method to 
compensate for the difference, as demonstrated success- 
fully by the Whipple group. We have tried to apply this 
method to some CANGAROO-I data sets, but did not 
find it to be effective. This is due to the fact that ADC 
signals have already been padded by electronic noise, 
which were generated only when high voltages were sup- 
plied to the PMTs, just like the original "hardware" 
adding utilized in r nost first-gene ration ACT systeni s 
Fruin fc JelievlfT968t IWekes et aLlHOTl I CawlevI 119931) . 
The other sources of noise described in the following sec- 
tions have larger effects and parameter distributions of 
background events match reasonably well between ON- 
and OFF-source regions after considering them. There- 
fore, software padding is not used in this paper. 

In 1997, we found that the discriminators which pre- 
cede the TDCs and scalers generate noise when they are 
fired, significantly affecting the ADCs on the same cir- 
cuit board. This results in inherent offsets in individ- 
ual ADCs in such events. This "ADC offset" noise can 
however be measured by using the ADC itself and by set- 
ting discriminator thresholds to be very low compared to 
their normal values (typically three photoelectrons). In 
the most affected channels, the offsets amount to about 
40 ADC counts, which correspond to about eight pho- 
toelectrons and thus distorts dimmer Cherenkov images 
significantly. We use the following formula to compen- 
sate for the ADC offset Oi instead of Equation ([1]): 



flj - Pi 



9i 



(2) 



From the pixels remaining after the above selections, 
isolated pixels, which are not adjacent to any other pixels 
containing signals, are further removed to extract clus- 
ters in the image. This process was traditionally designed 
to remove effects of the NSB noise. The remaining pix- 
els are used to define the image parameters described in 
Section [ 



4.3. Rejection of Electronic Noise Events 

The imaging camera consists of box units in which 
eight adjacent PMTs selected to have similar gains share 



the same high voltage. In early data, in particular, noise 
events in which all eight PMTs in a box unit were simul- 
taneously hit were frequently generated by the telescope 
drive system. These noise events, referred to as "box 
noise" , must be rejected as much as possible because they 
are compact and mimic gamma-ray events. We define a 
parameter h to reject box noise as follows: 

max rii 

b^- , (3) 



where rii is the number of TDC hits in the ith box of 
eight PMTs and n is the total number of TDC hits in an 
image. The box noise events have large values of b close 
to (or equal to) its maximum value of 1. In the standard 
analysis here, events of b greater than 0.8 are rejected. 
This rejection is very effective unless two or more boxes 
are triggered in this way simultaneously, although small 
Cherenkov images are also somewhat reduced 
rifice. About 15% of reconstructed events are rejected 
with this cut on average, but the fraction can be larger 
depending on the noise level. 

Using signals of the selected pixels, the image size (sum 
of the ADC values), the number of selected pixels (re- 
ferred to as A^hit though it is not the true number of 
TDC hits) , and the image centroid (first moment of the 
image) are calculated. Events of small size or A^^it are 
more affected by noise, and therefore events of size < 200 
ADC counts or Ni^n < 5 are rejected before the follow- 
ing imaging analysis. Also, images of the centroids more 
than 1?05 distant from the center of the camera are elim- 
inated. This is called the "edge cut" since part of such 
an image has possibly been lost out of the camera edge, 
resulting in the image shape being distorted. On av- 
erage, the size, A^hit, and edge cuts individually reject 
about 40%, 25%, and 55% of reconstructed events, re- 
spectively, and about 75% of reconstructed events are 
rejected by the three cuts applied all together. 

4.4. Imaging Analysis 

The atmospheric Cher enkov i maging a nalysis is based 
on the Hillas pa rameters ()Hilla£lfT 985: Wee kes et alll 19891 : 
iRevnolds et"all 'l993'l. which are the second moments of 
light intensity both parallel (length) and perpendicular 
(width) to the major axis of the image, the light concen- 
tration of the two brightest pixels (cone), the distance 
of the image centroid from the source (dis), and the ori- 
entation angle with respect to the major axis (alpha). 
They are also referred to as shape (width, length, and 
cone), location (dis), and orientation (alpha) parame- 
ters. Images of gamma-ray-induced air showers tend to 
have smaller width, length, and alpha, and larger cone 
than background hadronic showers, reflecting the differ- 
ent profiles of shower development and different source 
distributions in the field of view. In the case of a point 
source, the signal-to-noise ratio is improved by more than 
an order of magnitude by applying selections for gamma- 
ray-like events to the parameter values. 

The standard gamma-ray selection criteria used in this 
paper are simple and represented as follows: 



dr,l\u 



width 
length 
cone 
dis 
alpha 



^max 



(4) 
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Fig. 2. — Correlations between the best gamma-ray selection 
criteria giving the maximum q factors in the case of the Vela 1993- 
1995 analysis. Bottom left: Zmax (upper limit for length) vs. uimax 
(upper limit for width), top left: Cmin (lower limit for cone) vs. 
Wmax, and bottom right: Cn^in vs. imax- The star marks indicate 
the best combinations of the selection criteria and the lines indicate 
their Icr error regions. The different colors indicate different offset 
angles of the source from the center of the camera with the 0?1 
interval. 

The best combination of the upper and lower hmits in the 
above selection criteria is determined so as to maximize 
the quality factor q: 



q 



(5) 



width and length. The positions of the best shape selec- 
tion levels indicated by the star marks seem quite widely 
spread within the la error regions, especially for cone, 
since g is a discrete function of the selection levels by 
definition. 

The best shape selection levels for individual data sets 
are listed in Tabled For each data set, 5 x 10^ gamma- 
ray showers have been simulated, randomly injected into 
a circular area of 300 m radius around the telescope. 
Gamma-ray energies are also randomly selected in the 
range between 500 GeV and 100 TeV assuming that they 
follow a power-law differential spectrum or that with an 
exponential cutoff. The assumed spectral index and cut- 
off energy are listed in Table[3]together with the other in- 
put conditions in the simulations described in Section [2l 
The spectral index and cutoff energy for Vela are now 
known and taken from the reference given in the table, 
but the spectral indices for PSR B1706-44 and SN 1006 
are assumed to be the Crab-like value of —2.5. The zenith 
angles of gamma-ray injections are selected near the cul- 
minations of the objects at which most events were ob- 
served. In spite of the various differences in the obser- 
vation conditions, the best selection levels are consistent 
within the errors, of which the definition is the same as 
in Figure m except for the effect of the source offset from 
the camera center. 

The optimization of selection levels on location and ori- 
entation parameters is complicated if the effect of source 
extent is considered. Instead, we here take just simple 
criteria: the dis limits are common as dmin = 0?7 and 
c'lnax = 1°2 since their dependence on the source extent 
is not so strong, and amax = 10° for sources whose extent 
is comparable to the PSF, and for sources more extended 
the same ranges as used in the previous CANGAROO-I 
paper (Tanimori et al. 1998; alpha < 15°) are adopted. 



where e-y and ebg are the efficiencies of the selection cri- 
teria for gamma-ray and background events, respectively. 
Gamma-ray event samples are generated using Monte 
Carlo simulations, and for background event samples ac- 
tual OFF-source data are used. Only shape parameters, 
width, length, and eonc, are used in the optimization be- 
cause the distributions of location and orientation pa- 
rameters depend on the source distribution that is not 
known a priori. Nevertheless, the best selection criteria 
on the shape parameters depend on the source position 
in the camera, especially on the offset angle of the source 
from the camera center, which can change night by night 
owing to the shift of the tracking center. The average off- 
set angle is calculated in each data set and incorporated 
in the Monte Carlo simulations. 

Figure [2] shows correlations between the best gamma- 
ray selection levels on the shape parameters with the 
OFF-source data of the Vela 1993, 1994, and 1995 ob- 
servations used as the background samples. The best 



combination of w 



max; •'max; 



and Cmin for each offset angle 



of a point source from the camera center is indicated by a 
star mark, and the closed line around it represents the la 
error region, which is defined so that q > qmax — crq^._^^. It 
is clearly seen that the more distant the source position 
is from the camera center, the larger the upper limit to 
length /max is, but Wmax a-nd Cmin are less dependent on 
the source offset. The relative errors of Cmin are larger 



4.5. Estimation of Integral Gamma-Ray Fluxes 

The integral gamma-ray fluxes obtained in this paper 
have been estimated using the following formula: 



F{> Sth) 



expected 



-Fo{> Etk) 



(6) 



where i?th is the threshold energy for gamma rays, which 
is defined as the modal energy of selected gamma rays, 
-^excess is the number of excess events, Fq is the assumed 
integral gamma-ray flux from the source, and iVoxpcctod 
is the expected number of gamma rays represented as 
follows: 



expected 



Aoet 



k{E)dE, 



(7) 



where /o = —dF^/dE is the assumed differential gamma- 
ray flux from the source, -Emin and Emax are the mini- 
mum and maximum energies of gamma rays generated in 
the simulations, respectively, Aq = Trr^j^^ cos^ is the area 
within which simulated gamma rays were injected, with 
the maximum horizontal distance between the telescope 
and the shower axis rmax = 300 m and the zenith angle 

^/A^inputi 



9, e is the acceptance defined as e = A^si 



in which Mnput and N, 



surviving 



are the numbers of input 
gamma rays and gamma rays surviving after the selec- 
tions, respectively, and t is the observation time. 
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TABLE 2 

Best gamma-ray selection criteria for each data set. 



Target 


Year 


Offset^ 




9m ax 




^max 


^max 


^min 


PSR B1706-44 


1993 


0?1 


1 


549 ± 


012 


0?0883+ o""^ 

-or 0036 


0? 288+0^0" 


0.17t0;18 




1994 


0?1 


1 


567 ± 


015 


onr>o^o-t-0°0057 

0?0873^ o 

-or 0033 


0?288+0o025 
-0?008 


„ „ . -1-0 03 




1995 


0?1 


1 


563 ± 


014 


0?0896^ o 

-or 0049 


0?292+"o 

-0?012 


„ ^ --1-0 20 

O-lSIo.Is 




1997 


0?1 


1 


748 ± 


016 


0?0856^ o 

-or 0058 


„^„„ .i+0°017 
0?294+ o" 
-oroii 


0-08Io.08 




1998 


0?1 


1 


995 ± 


025 


^^n^o-. ^^-t-0°0047 

0?0810+OoOO^J 


.-.0^ + 0° 015 

0-300^0? 026 


^ -. ^-1-0 25 
0.12tO-25 


Vela 


1993, 1994, 1995 


0?2 


1 


560 ±0 


007 


0?0886«i°o- 


00296+o?oi2 








0?4 


1 


466 ±0 


006 


0°OQll+°°°°^'^ 
'J'jy^^_0?0040 


nooon + O. 016 


n 20+013 




1997 


0?1 


1 


635 ±0 


008 


0? 0929+0- 

-0?0066 


0? 296+0^016 

-0?009 


'^■^■-'-0.23 






0?4 


1 


482 ±0 


007 


0?0953+°-°036 

-0?0058 


0? 331+0-028 

-0?014 


in+0-22 


SN 1006 


1996 


0?3 


1 


557 ±0 


022 


0?0962+°-°°^2 

-0?0033 


0? 309+0-025 

-0?012 


r, gq + 0.08 

*J-^'^-0.23 






0?7 


1 


413 ±0 


018 


0? 0957+0- 

-0?0041 


0? 393+0-028 

-0?065 


1^+013 

"■^'^-a.is 




1997 


0?3 


1 


555 ±0 


010 


0? 0924+0- 0030 

-0?0063 


0?321+0-0i4 

-0?021 


n+o-21 

'J-l'^-0.13 



Offset angle of the source from the center of the camera. 



TABLE 3 

Input conditions in the Monte Carlo simulations for each data set. 







Cutoff Energy 


Zenith Angle 


Reflectivity 


Target 


Year 


Spectral Index (TeV) 


(deg) 


Light Guide (%) Number of PMTs 



PSR B 1706-44 


1993 


-2.5 




15 


no 


45 


224 




1994 


-2.5 




15 


no 


45 


224 




1995 


-2.5 




15 


no 


45 


256 




1997 


-2.5 




15 


no 


75 


256 




1998 


-2.5 




15 


no 


55 


240 


Vela 


1993, 1994, 1995 


-1.45=' 


13.8'' 


15 


no 


45 


224 




1997 


-1.45=" 


13.8'' 


15 


no 


75 


256 


SN 1006 


1996 


-2.5 




10 


yes 


45 


256 




1997 


-2.5 




10 


no 


75 


256 



^ lAharonian et al.l ll2006al). 



5. RESULTS 

The numbers of events remaining after the selections 
described above are summarized in Table |4] for all data 
sets used in this paper. The selections listed in the ta- 
ble header have cumulatively been applied to each data 
set from left to right. The statistical significances of the 
excess events obtained after applying all select ions have 
been calculated using the formula obtained bv lLi fc Mai 
{1983) utilizing the likelihood ratio. Note that after 
the noise rejections the numbers of ON- and OFF-source 
events match "with each other -within a fe^w percent ex- 
cept for the 1998 PSR B1706-44 data. Details are given 
in the follo"wing subsections. 

5.1. PSR B1706-U 

The five data sets of PSR B 1706-44 have been reana- 
lyzed separately since the 1993 data used in the previous 
CANGAROO-I paper should independently be analyzed 
for comparison, and the mirror reflectivities and camera 
configurations of the other data sets listed in Table [1] all 
differ. The average offset of the object from the cam- 
era center is about 0?1, -which has been assumed to be 
common for all data sets. The gamma-ray selection crite- 



ria listed in Table [2] are optimized for this common offset, 
but the selection levels are different from data set to data 
set since the OFF-source data of each data set have been 
used as the background samples in each optimization. 

The alpha distributions of ON- and OFF-source events 
remaining after the standard selections have been plotted 
for individual data sets in Figure [31 No statistically sig- 
nificant excess of ON-source events has been found near 
alpha = 0° in any data set, and the significant excess 
from the 1993 data previously reported has not been 
reproduced in this analysis. The ON- and OFF-source 
alpha distributions of each data set agree "well except 
for the 1998 data, in -which the ON-source distribution is 
some-what higher than the OFF-source distribution in the 
"whole range. A possible reason for this is a temporary 
decrease of the mirror reflectivity due to de"w condensa- 
tion since it "was humid in this period and more OFF- 
source data "was taken in the hours before dawn. This 
diffuse excess is not gamma-ray-like considering the dis- 
tributions of the shape parameters of the excess events, 
and a conservative upper limit is calculated here leaving 
the background mismatch as it is. Setting Cmax — 10° 
for a relatively point-like source and assuming a Crab- 
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TABLE 4 

Summary of numbers of events remaining after various selections. 



Target 


Year 


ON / OFF 


Raw^ 


Noise 


Shape° 


T 4-- d 

Location 


Orientation" (amax) 


Excess 


Significance 


PSR B1706-44 


1993 


ON 


180105 


57587 


14942 


8669 


1090 (10°) 


-67 


-1.4(T 






OFF 


167397 


56085 


14570 


8469 


1157 (10°) 








1994 


ON 


62616 


20602 


5014 


2763 


366 (10°) 


-24 


-0.9(T 






OFF 


63818 


21186 


5094 


2863 


390 (10°) 








1995 


ON 


159487 


19865 


4454 


2445 


299 (10°) 


-31 


-1.2ct 






OFF 


138384 


20508 


4775 


2664 


330 (10°) 








1997 


ON 


114862 


22541 


3736 


1994 


230 (10 ) 


— 2 


— O.lcr 








112291 


22651 


Q7Q9 


2062 










1998 


ON 


91547 


18662 


2422 


1280 


160 (10°) 
138 (10°) 


22 


1.3(7 






OFF 


84319 


16972 


2082 


1090 






Vela 


1993, 1994, 1995 


ON^ 


420310 


119974 


29710 


16970 


2343 (10°) 


298 


4.5(7 






OFF^ 


400418 


120206 


28867 


16512 


2045 (10°) 










ons 


420310 


119974 


35190 




10122j 


501 


3.6(7 






offs 


400418 


120206 


34243 




962lj 








1997 


ON^ 


142158 


36064 


8039 


4421 


589 (10°) 


84 


2.5(7 






OFF^ 


132171 


35292 


8004 


4309 


505 (10°) 










ons 


142158 


36064 


9862 




2624j 


75 


1.0(7 






OFPS 


132171 


35292 


9772 




2549j 






SN 1006 


1996 


onI^ 


33524 


2105 


564 


216 


32 (15°) 


-5 


-0.6(7 






OFF*^ 


27171 


2122 


581 


213 


37 (15°) 










ON' 


33524 


2105 


696 


274 


46 (15°) 


2 


0.2(7 






OFF' 


27171 


2122 


728 


266 


44 (15°) 








1997 


on'' 


231398 


57143 


13572 


6915 


1246 (15°) 


-94 


-1.8a 






off'' 


226916 


57112 


14016 


7199 


1340 (15°) 







^ Number of reconstructed events (at least two adjacent pixels containing signals). Number of events remaining after the noise rejection described 
in Section 14.31 and the discharge noise rejection described in Section [5T2I for the SN 1006 1996 data.° Number of events remaining after the shape 
selections based on width, length, and conc.^ Number of events remaining after the location selection based on cZis.*^ Number of events remaining after 
the orientation selection based on alpha.^ Analyzed with respect to the CANGAROO-I position (a — 8*^35^"42^, S — —45° 17' (J2000)).^ Analyzed 
with respect to HESS J0835-455 (a ^ S^SS'^OO^, S ^ -45°36' (J2000)).*^ Analyzed with respect to the NE rim with the selection criteria optimized 
for the source offset by 0?3 from the camera center.^ Analyzed with respect to the NE rim with the selection criteria optimize d for the source offset 
by 0?7 from the camera center.-' Selected so that the arrival direction obtained using the source PDF method llYoshikoshilll996l1 is within the angular 
distance of 0?8 from HESS J0835-455. 



like spectral index of —2.5 in the simulations, the upper 
limits to the integral fluxes at the 95% confidence level 
(CL) have been estimated using the method des cribed 
in Section [45] and the method of iProtherol (| 19 841 ). The 
results are as follows: 

F93(> 3.2±1.6TeV) < 8.03 x 10"^^ photons cm^^ s~\ 
^93-94(> 3.2±1.6TeV) < 6.06 x 10"" photons cm^^g-i^ 
F95(> 3.2±1.6TeV) < 8.85 x 10"" photons cm"^ s"\ 
F97(> 1.8±0.9TeV) < 4.08 x 10"" photons cm"^ s"\ 

F98(> 2.7±1.4TeV) < 1.29 x 10"" photons cm"^ s"\ 

where -F93-94 is the integral flux calculated adding the 
1993 and 1994 data sets in which the observation condi- 
tions are the same. The different threshold energies for 
the above upper limits are mostly due to the different 
mirror reflectivities listed in Table [3l The errors of the 
threshold energies are systemat ic errors, define d in the 
same way as described by Murai shi et al.l ()200Cl[ ). In the 
calculation of the flux upper limit of the 1997 data, only 
7.9 hr of the ON- and OFF-source data, which corresponds 
to 43% of the data used elsewhere, have been used since 
blue optical filters were tested in the rest of the observa- 
tions and their effect has not accurately been determined. 
The upper limits are plotted in Figure[4l They are higher 
than, and therefore not inconsistent with , the H .E.S.S. 
upper Hmits reported by I Aharonian et all (|2005a[ ). 



5.2. SN 1006 

The SN 1006 data taken in 1996 and 1997 have 
been analyzed with respect to the northeast (NE) rim 
(a = 15''03°'5y S = -41°45'30" (J2000)), from which 
iTanimori et al.l (|1998f ) have claimed the detection of TeV 
gamma rays. The observations were made with aver- 
age offsets of 0?3 (June 1996 and 1997) and 0?7 (1996 
April) of the NE rim from the camera center in order 
to move the bright star f3 Lup (magnitude 2.68) out of 
the field of view of the camera. Because of these rela- 
tively large offsets, the best gamma-ray selection criteria 
for SN 1006 listed in Table [2] are looser than the others 
following the tendency indicated in Figure [H In the case 
of the 1996 data, the looser criteria are also caused by 
the light guides^^ which were used only in this period. 
In the case of CANGAROO-I, the collection efficiency of 
the light guides is thought to be worse than originally ex- 
pected owing to difficulties in placing them in proximity 
to the photocathodes which resulted in discharge noise 
generated between the light guides and photocathodes, 
and the Cherenkov event rate in fact decreased after their 
installation. Here, to generate simulation events, we as- 
sume the light guide efficiency to be 45%, which is almost 
the same as the filling factor of the photocathodes. As in 
the case of PSR B1706-44, the Crab-like spectral index 
—2.5 is assumed in the simulations. 

Light guides are utilized to detect Cherenkov photons other- 
wise incident on the gaps between photocathodes of PMTs but at 
the same time serve to somewhat smear Cherenkov images. 
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Fig. 3. — Distributions of the orientation parameter alpha for 
tile PSR B1706-44 data sets. From tlie top: distributions of the 
1993, 1994, 1995, 1997, and 1998 data. The histograms with error 
bars are the ON-source distributions, and the gray solid histograms 
are the OFF-source distributions. 



As mentioned above, the 1996 data were contaminated 
by discharge noise, in which the hardware trigger rate 
suddenly increased, especially in humid conditions, with 
only one or two PMTs having very large ADC values 
in one image. Such images are compact and have large 
sizes, and therefore, mimic gamma-ray images. We here 
take the following simple method to reject the discharge 
noise events. 




10 

Energy (TeV) 

Fig. 4. — Integral gamma-ray fluxes from PSR B1706-44 ob- 
tained using the CANGAROO-I data. The 95% CL upper limits 
to the integral flux from PSR B1706-44 obtained from the 1993, 
1993 -I- 1994, 1997, and 1998 data are indicated by the filled circle, 
filled square, filled triangle, and filled upside-down triangle, re- 
spectively. The ope n circle indicates the integral flux estimated by 
IKifune et al.l II1995I) using the 1993 data. The horizontal bars rep- 
resent the systematic errors of the thresho ld energies. The 99% CL 
upper limit curves obtained by H.E.S.S. HAharonian et al.|[2005al) 
are indicated by the thin solid, dashed, and dotted lines. The thick 
solid line is the integral spectrum converted from the preliminary 
CANGAROO-II result IIKushida et al.ll200^). The o pen square 
and flUed sta r indic ate the results from lChadwick et ah ( 1998) and 
IRowell et al.l II1998I ). respectively. 

1. For each pixel, the rate in which the pixel has the 
maximum ADC value in the image is calculated. 

2. Pixels with rates higher than 0.01 Hz are marked 
as noisy pixels night by night (the rate of normal 
pixels is about 0.001 Hz). 

3. Discharge noise events are identified and rejected 
if the pixel having the maximum ADC value in the 
image coincides with one of the marked pixels for 
the night. 

This discharge noise cut is more robust than the old 
method described in Section 16.21 and more than 75% of 
shower events remain after this cut. The numbers of ON- 
and OFF-source events after the cut match quite well and 
the difference between them is only 0.8%. 

The 1997 data are free from the discharge noise since 
the light guides were removed in 1996 November. How- 
ever, the number of ON-source events after the other 
noise cuts described in Section 14.31 is more than a few 
percent larger than that of OFF-source events. This mis- 
match is probably due to the difference of brightness in 
the field of view between ON- and OFF-sources since an- 
other bright star, k Cen (magnitude 3.13), was inevitably 
included in the ON-source field of view. To reduce the ef- 
fect of bright stars, pixels having scaler values greater 
than 100 counts ms^^ were removed from the image in 
the analysis of the 1997 data. As a result, the numbers 
of ON- and OFF-source events match quite well as shown 
in Table H 
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Fig. 5. — Distributions of the orientation parameter alpha for 
the SN 1006 data sets. The alpha values have been calculated 
with respect to the NE rim of the SNR. Top: distributions of the 
1996 data with the best selection criteria optimized for the 0?3 
offset from the camera center, middle: distributions of the 1996 
data with the best selection criteria for the 0?7 offset, and bottom: 
distributions of the 1997 data with the best selection criteria for 
the 0?3 offset. The histograms with the error bars and the gray 
solid histograms are ON- and OFF-source distributions, respectively. 



Figure [5] shows the alpha distributions of the SN 1006 
data sets. No statistically significant gamma-ray signal 
has been found around alpha = 0° from the 1996 and 
1997 data for both 0?3 and 0?7 offsets from the cam- 
era center. Therefore, the previous CANGAROO-I re- 
sults have not been reproduced in this analysis. Set- 
ting Omax = 15° , whic h is the same definition as used by 
iTanimori efall (|1998f ). the 95% CL upper limits to the 
integral flux have been calculated as follows: 

^96 0?3(> 3.0±1.5TeV) < 1.20x 10"^%hotons cm^'^s'^, 

FgggO.^(> S.Oil.STcV) < 2 . 4 4 X l Q ^ ^ ^ p ho t o US cm^^s"^ 

^97,o?3(> l-8±0.9TeV) < 1.96x IQ-^^ photons cm^^g-i^ 

The difference between the threshold energies is mostly 
due to the different mirror reflectivities listed in Table [3l 
The errors of the threshold energies are systematic er- 
rors, as defined by Muraish i et al.l (I200GD . These upper 
limits are plotted in Figure [S] They are higher than and 
not inconsiste nt with the H.E.S.S. up per limits on the 
SN 1006 flux (jAharonian et al.ll20051 
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Fig. 6. — Integral gamma-ray fluxes from the NE rim of SN 1006 
obtained using the CANGAROO-I data. The 95% CL upper lim- 
its to the integral fluxes are indicated by the filled circle (1996 
data with the selection criteria optimized for the 0?3 offset) and 
the filled square (1997 data). The open circle and open square 
are the integral fluxes estimated by Tanimori et al. (1998) using 
the 1996 and 1997 data, respectively. The thick solid line is the 
integral sp ectrum converted fro m the CANGAROO-I differential 
spectrum I ITanimori et al.l 120011) . The horizontal bars represent 
the systematic errors of the threshold en ergies. The 99 .9% CL 
upper limits obtained by H.E.S.S. (Aharoni an et aLll2005tD are in- 
dicated by the gray band for a range of assumed photon indices 
(2—3). T he open triangle i ndicates the preliminary CANGAROO- 
II result IIHara et aLlljoOlT) . 



5.3. Vela 

The analysis procedure has been applied to the Vela 
data, first with respect to the CANGAROO-I posi- 
tion {a = 8'^35'"42", S = -45°17' (J2000)), which 
is 0?13 offset from th e Vela pulsar to the southeast 
(|Yoshikoshi et al.lll997D . The 1993-1995 data and 1997 
data have been analyzed separately since the mirror re- 
flectivity was improved in 1996. The average offsets of 
the object from the camera center are about 0?2 (1993- 
1995) and 0?1 (1997), for which the shape selection cri- 
teria listed in Table [5] have been optimized. The ob- 
tained alpha distributions are shown in Figure [Tj Setting 
Omax = 10°, the statistical significances of the 1993- 
1995 data and 1997 data are 4.5a and 2.5cr, respectively. 
The former is smaller than the previous result of 5.2a 
(jYoshikoshi et al. 1997), but the signal is still at a sig- 
nificant level. Note, however, that the CANGAROO-I 
position was found in the previous analysis after search- 
ing for the h ighest peak in th e 2° x 2° area around the 
Vela pulsar (jYoshikoshi et al.]ll997t l. and therefore the 
above values are "pre-trial" statistical significances. 

Considering the 0?18 PSF (half width at half max- 
imum) of CANGAROO-I ( Yoshikoshi et al.lll99l . the 
CANGAROO-I position does not coincide with the peak 
position of HESS J0835-455 detected b y H.E.S.S., which 
is sho wn at the top left of Figure [5] (jAharonian et al.l 
'2006a). However, the acceptance e is not uniform over 
the field of view and the morphology could significantly 
be distorted by this non-uniformity, especially for the 
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Fig. 7. — Distributions of the orientation parameter alpha for the Vela data. The alpha values have been calculated with respect to the 
CANGAROO-I position that is offset by 0?13 from the Vela pulsar to the southeast. The left and right figures have been obtained from 
the 1993-1995 data and the 1997 data, respectively. The ON- and OFF-source distributions are represented by the histograms with the error 
bars and the solid histograms, respectively. The histograms in the insets are the distributions of excess counts. 



CANGAROO-I imaging camera which had a relatively 
small field of view, about 3° across. The top right of 
Figure [S] is the averaged acceptance distribution in the 
Vela pulsar region for the 1993-1995 observation pe- 
riod, obtained considering the night-by-night shifts of the 
tracking center in the camera and rotation of the cam- 
era with respect to the sky view due to the alt-azimuth 
mount. The expected morphology for CANGAROO-I is 
obtained from the H.E.S.S. morphology reweighted ac- 
cording to the CANGAROO-I acceptance and is shown 
at the bottom left of Figure [51 where the accept ance dis- 
tribution of H.E.S.S. (Aharonian et ,al.| |2006b() has not 
been considered but is almost uniform within the field of 
view compared to the CANGAROO-I acceptance. In this 
acceptance-reweighted morphology, the head of the "sea- 
horse" shape, which corresponds to the CANGAROO-I 
source as shown in the bottom right of Figure [H is more 
emphasized than in the original one. The CANGAROO-I 
morphology is obtained using the source proba bility den- 
sity function (PDF) method (|Yoshikoshilll996( l. in which 
a PDF for the source position with respe ct to the im- 
age centroid, ma jor axis, and asymmetry ()Punchlll993t 
lYoshikoshil 119961 ) is defined using Monte Carlo simula- 
tions and a source distribution is made by adding up all 
PDFs of gamma-ray-like events. The angular distance 
between the peak of the acceptance-reweighted H.E.S.S. 
morphology and the CANGAROO-I position is 0?05, 
which is comparable to the la error of the CANGAROO- 
I position - 0?04 (|Yoshikoshil |1996[ ). Therefore, the 
TeV gamma-ray signal detected by CANGAROO-I is 
possibly part of HESS J0835-455. However, the above 
acceptance-reweighted morphology is still more extended 
than the CANGAROO-I source. This is possibly due to 
the energy dependence of the gamma-ray morphology: 
the threshold energy of the H.E.S.S. observations of Vela 
is 450 GeV which is an order of magnitude lower than 



that of CANGAROO-I (see Section [O] for further dis- 
cussion) . 

The integral flux from HESS J0835-455 has been esti- 
mated using the method described in Section 14.51 This 
H.E.S.S. source is diffuse extending up to the radius of 
0?8 centered on position II (a ^ S^SS ^OO", S = -45°36' 
(J2000)) defined by Aharo nian et all ()2006aD . The dis 
and alpha cuts are not useful to estimate the flux of the 
same circular area as used by H.E.S.S., and instead we 
must estimate the arrival direction of each event as in 
the case of stereoscopic observations. Here, we estimate 
the arrival direction as the most probable point of the 
source PDF described above. The best gamma-ray selec- 
tion criteria have been obtained for the position offset by 
0?4 from the camera center since the position of the max- 
imum emission of HESS J0835-455 is located around it 
on average. Simulation events have been generated with 
a power-law spectrum of th e index of —1.45 and an ex- 
ponential cutoff at 13.8 TeV (|Ah aronian et al.ll2006a[ ). In 
the case of the 1993-1995 data, the numbers of ON- and 
OFF-source gamma-ray-like events falling into the circu- 
lar area are 10,122 and 9621, respectively. The excess 
of 501 events corresponds to 3.6a, which is less signifi- 
cant than that for the CANGAROO-I position, owing to 
the inclusion of more background events in the extended 
area. The acceptance e depends on the source distri- 
bution in the field of view, but first we simply assume 
that the source is uniformly distributed in the circular 
area and e is averaged over it. The integral flux from 
HESS J0835-455 (0?8 radius centered on position II) has 
thus been obtained as follows: 

-F93-95(> 4.0 ± 1.6 TeV) = (3.28 ± 0.92) x lO"^^ 

photons cm^^ s^^. 

The stated error in the t hreshold energy i s the system- 
atic error, as defined by iMuraishi et al.l (|2000f ) except 
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Fig. 8. — Sky maps of the Vela p ulsar region. Top left: distribution of excess counts obtained by H.E.S.S. with the threshold energy 
of 450 GeV IIAharonian et aT1l2006an . top right: averaged acceptance distribution in the CANGAROO-I observations from 1993 to 1995 
normalized so that the maximum value is e qual to 1, bottom left: H .E.S.S. excess distribution rowoightod according to the CANGAROO-I 
acceptance, where the H.E.S.S. acceptance (Aharonian et al." 2006b') has not been considered but is almost uniform within the field of view 
compared to the CANGAROO-I acceptance, bottom right: distribution of excess counts obtained by CANGAROO-I with the threshold 
energy of 4TeV using the source PDF method (iYoshikoshiiil996il . in which the PSF size is 0?18 (half width at half maximum). The rainbow 
scales at the right are in counts per 0?02 X 0?02 area for the H.E.S.S. distributions and in counts deg~^ for the CANGAROO distribution. 
The star mark at the center of each sky map indicates the position of the Vela pulsar. The filled circle indicates the peak position of the 
CANGAROO-I excess, whereas the filled and open squares indicate the position of the center of gravity of the H.E.S.S. excess and the peak 
position of the acceptance-reweighted H.E.S.S. excess, respectively. The small open circle around the CANGAROO-I position represents 
the 2(T location e rror (la ~ 0?04 ). The big dashed circle of the 0?8 radius indicates the integration region for HESS J0835— 455 defined by 
H.E.S.S. (.Aharonian et al.|[2006al) . 



for the uncertainty of the spectral index, whereas the er- 
ror of the flux is only statistical. The threshold energy of 
4.0 TeV is higher than t he previous estimation of 2.5 TeV 
(jYoshikoshi et al.lll997t l. owing to the harder spectrum 
used in the simulations this time. For HESS J0835-455, 
the two-dimensional Gaussian profile of the diffuse emis- 
sion has been given by H.E.S.S. The average acceptance 
can more realistically be calculated using this profile, as 
weights and the integral flux obtained using this accep- 



tance is 

^93-95,g(> 4.0 ± 1.6 TeV) = (3.05 ±0.86) x 10"^^ 

photons cm^^ . 

The difference between the above two fluxes is small and 
well within the statistical errors. 

The same calculation has been done also for the 1997 
data: 2624 ON-source and 2549 OFF-source events re- 
mained after the selections and the excess of 75 events 
corresponds to l.Ocr. Since this excess is not statistically 
significant, the 95% CL upper limits to the integral flux 
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2. The gamma-ray selection criteria were defined as 
follows: 
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Fig. 9. — Integral gamma-ray fluxes from the Vela pulsar re- 
gion (0?8 radius centered on HESS J0853-455) obtained from the 
CANGAROO-I data. The open circle represents the integral flux 
from the 1993 to 1995 data, and the flUed circle represents the 95% 
CL upper limit to the flux from the 1997 data. The horizontal bars 
represent the systematic errors of the threshold energies. The solid 
line is the integral spectrum converted from the differential spec- 
trum of HESS J0853-455 measured b y H.E.S.S. assumi ng a power 
law with an exponential cutoff ( Aharonian et ani2006al) . 

have been calculated as follows: 

Fg7(> 2.5 ± 1.0 TeV) < 8.87 x 10"^^ photons cyr-'^s-'^, 

Fq7,g{> 2.5±1.0TeV) < 8.33 x 10"^^ photons cm"2s-i_ 

These fluxes are plotted in Figure [5] as well as the inte- 
grated H.E.S.S. spectrum fjAharonian ct al..2006_a). The 
estimated fluxes are all consistent with the H.E.S.S. spec- 
trum within statistical errors. 

6. DISCUSSION 

The results of the reanalyses described in the previ- 
ous section are consistent with the recent H.E.S.S. re- 
sults except for the morphological difference in the Vela 
pulsar region, which is possibly due to the difference 
of the threshold energy between the two experiments. 
In other words, the previous CANGAROO-I results for 
PSR B1706-44 and SN 1006 claiming the significant de- 
tections of TeV gamma rays have not been reproduced 
in the reanalyses. Therefore, we have tried to apply the 
same analysis methods as used before to the PSR B1706- 
44 and SN 1006 data to determine why the old gamma- 
ray signals were obtained, although not all details of the 
original analysis are still available. 

6.1. Reproducibility of the Old PSR B 1706-44 Results 

Details of the analysis used to obtain the old 
PSR B 1706-44 res ults are not given in the paper by 
iKifune et al.l (fl995h . The refore, we have used the anal- 
ysis method described bv iTamura et al.l (1994) instead. 
The known differences of the method from the standard 
analysis described in Section |4] are listed below. 

1. Some noise cut levels are different: size > 100, 
A^hit > 4, and b (box noise parameter) < 0.7. 



0?045^ 



r\o A size 
"■^1500 



0?01 < width < 0?14 

0?1 < length < 0?6 

length < 0?2-|- 

0?2 < dis < 0?96 
length < dis 

alpha < 6? 75. 



In the old analysis, the above conditions were applied 
only to the PSR B1706-44 data taken in 1993 August, 
and the details used in the analysis of the other 1993 
data are not given. Therefore, we applied the above con- 
ditions only to the 1993 August data, but did not find 
any significant alpha peak around 0°. 

Subsequent examination of old internal collaboration 
documents revealed the following facts about the 1993 
August analysis. 

1. Only seven of the nine ON-source runs taken in 1993 
August were used in the analysis. 

2. Calibrated positions of the tracking center (Sec- 
tion l4.ip were not used since the calibration method 
had not yet been established at that time. Instead, 
the position of the tracking center was fixed in the 
camera at the position of (x, y) — (0?1, — 0?125), 
which is offset from the calibrated tracking centers 
by 0?13 on average. 

3. In the image cleaning, the threshold for ADC values 
was fixed at 5 ADC counts. 

4. Noise cut levels different from those in the standard 
analysis were used: 100 < size < 5000, TVhit > 4, 
and b < 0.7. Also, the effect of the ADC offset 
noise (see Section l4?2l) was not corrected. 

5. Another noise rejection called the "region cut" was 
applied. Events with centroids located within the 
areas of "noisy" box units were rejected in this cut. 
The "noisy" boxes were determined run by run and 
the gross number of the "noisy" boxes was eight in 
1993 August. 

6. The gamma-ray selection criteria were defined as 
follows: 



0?01 < 



0?1 < 



0?2 
length 



width 
width 
length 
length 

dis 

dis. 



0?25 
0?1 -f 
0?6 
0?2 + 
0?96 



no HQ size 



"■^"^1500 



The selection levels were modified slightly run by 
run, but no further information beyond this is avail- 
able. 

One can note that s ome of the above sele ction levels are 
different to those of iTamura et al.l (|1994l ). There seems 
to be sor ne small changes in the analysis betw e en th e 
papers of ITamura et al.r(|1994[ ) and IKifune et al.l 1199.5!). 

We have applied the above conditions to the 1993 Au- 
gust data and the obtained alpha distributions are shown 
in Figure [TOl Some excess of the ON-source events has 
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Fig. 10. — Distributions of the orientation parameter alpha 
for the PSR B1706-44 data taken in 1993 August. The analysis 
method is based on the descriptions in the old internal documents 
to try to reproduce the result obtained by Kifune et al. (1995). 
The histogram with the error bars and the gray solid histogram 
are the ON- and OFF-source distributions, respectively. The num- 
ber of events remaining af ter selections is similar to that in Figure 1 
(d) of lKifune etaLl llT995h . although the bin sizes of the histograms 
differ. 



been found near alpha = 0° and its statistical signifi- 
cance is 2.8(7 taking the events of alpha smaller than 10°. 
However, the excess is n ot as sharp as t he old result. Re- 
garding the significance, iKifune et al.1 ()1995[ ) used a dif- 
ferent definition, in which the background level was es- 
timated from the ON-source events of alpha > 30° since 
their alpha distributions appeared to be flat. Following 
this definition, the significance is enhanced by a factor of 
•\/2/(l + a) ^ 1.3, where a = 1/6 is the normalization 
factor . It is also found from Figure 1 oflKifu ne et al.l 
()1995( ) that the amount of the 1993 August data used in 
our reanalysis is roughly 1/5 of the data used in the old 
analysis. If we simply scale the total significance using 
the above factors, it turns out to be 2.8(7 x 1.3 x \/5 ^ 8a, 
which is not so different from the original claim of 12cr de- 
spite the very rough estimation. However, we think that 
this excess is not due to gamma rays from PSR B1706-44 
since the tracking center (0?1, — 0?125) in the camera is 
not very close to the calibrated tracking centers and the 
Cherenkov images have been left affected by the ADC 
offset noise. Only a broad alpha peak smeared by the 
source offsets and the ADC noise would be expected if it 
were really due to gamma rays. 

A possible reason for the apparently significant 
gamma-ray-like signal in the previous analysis is that an 
initial excess has been enhanced by the complicated se- 
lection levels, and/or their run-by-run modifications but 
the consequent statistical penalty for the number of de- 
grees of freedom used has been neglected. For example, 
the "region cut" was introduced in the name of noise 
rejection but the "noisy" boxes were not selected inde- 
pendently of the gamma-ray-like signal. If the 8 x 0.4 ~ 3 
"noisy" boxes are randomly selected among the total 



number of boxes, 28 boxes per camera x 7 runs x 0.4 
~ 78, the number of possible combinations in the "re- 
gion cut" is (J^) = 76, 076, where the factor 0.4 is the 
ratio of the effective camera area after the dis selection 
of dis < 0?96. The total gamma-ray-like signal (number 
of gamma-ray-like excess events) z of a combination is 
represented as 



78 



(8) 



where Xi E X = {xi, X2, xrs} is the signal of 

the ith box (number of excess events of which the cen- 
troids are located within the box area) and yj € X is 
also a signal of a box but the index is renumbered for re- 
jected boxes. Assuming that Xi is normally distributed 
with mean and standard deviation cr, any z is also nor- 
mally distributed with mean and standard deviation 
a' = — 3 = ^/75a. However, if the maximum z 

is chosen among all combinations, Zmax is no longer dis- 
tributed around 0, and in this case, the distribution has 
a positive offset of about la'. These kinds of offsets are 
accumulated by the number of degrees of freedom, and 
the signal can be made apparently significant just from 
the background fiuctuations. Therefore, we suspect that 
the very complicated conditions in the old PSR B1706~44 
analysis described above could supply sufficient degrees 
of freedom to generate the reported signal. 

We have also investigated how much the significance 
can be enhanced in the real data by varying the region 
cut. The selection criteria used in the old analysis have 
been adopted except the region cut, and the number of 
rejected boxes in varied region cuts has been fixed at 8. 
In the real application, the above approximation using 
the factor 0.4 is not possible because the region cut is 
discrete, and the total number of boxes increases to 147, 
which is derived from 21 boxes per camera fully or par- 
tially contained in the dis selection area, multiplied by 
the 7 runs. The dis selection is the same in both, but 
they differ in how to count boxes partially contained in 
the dis selection area. The number of boxes per camera 
before region cuts is 28 x 0.4 = 11.2 in the toy model, 
and 4 (fully contained) -1-17 (partially contained) =21 
in the real application. A total of 10^ combinations has 
randomly been sampled from the all (""^g^) ~ 4.5 x 10^^ 
combinations and their significances have been calculated 
after applying the corresponding region cuts. They are 
normally distributed around a mean significance of 2.4cr 
with a root mean square of 0.24(T. The maximum signif- 
icance is 3.5(7 with this number of samples. The signifi- 
cance obtained without using the region cut is 2.5(7 and 
thus the increase of the significance with the region cut 
actually used in the old analysis is 0.3tT of the total 2.8a, 
which is in fact a positive enhancement, although smaller 
than the possible maximum offset estimated above. How- 
ever, the probability of obtaining a significance greater 
than 2.8(7 with any region cut of 8 boxes is 6.4%, which 
is small to be interpreted by chance. 

6.2. Reproducibility of the Old SN 1006 Results 

Deta ils of the old SN 100 6 analysis are not given in 
full in iTanimori et al.l |1998) , an d so we have followed 
the analysis method described bv iKameil ()1998l) instead. 
The gamma-ray selection criteria were as follows: 
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1. for the 1996 April data, 

0?04 < width < 0?17 

0?10 < length < 0?48 
0.64 < cone 
0?1 < dis 
length < dis 

alpha < 15°, 

2. for the 1996 June data, 

0?04 < width < 0?17 

0?10 < length < 0?38 

0?4 < dis < 1?2 
3.3length~ 1?0 < dis 

alpha < 15°, 

3. for the 1997 data and events of size < 2000, 

0?04 < width < G?13 

-0?23 + 0?131ogsi2e < length < 0? 07 + 0? 084 log size 

0.23 log size < eonc < 0.42 + 0.14 log size 

0?5 < dis < 1?1 
length < dis 

alpha < 15°, 

4. for the 1997 data and events of size > 2000, 

0?04 < width < 0?13 
0?2 < length < 0? 07 + 0? 084 log size 
0.72 < cone < 0.42 + 0.14 log size 
0?5 < dis < 1?1 
length < dis 

alpha < 15°, 

where the defi nition of cone i s diffe rent from that orig- 
inally used by IWeekes et al.l (119891 ) and is the fraction 
of the light detected by the camera that is contained in 
the brightest half of the triggered pixels. The other con- 
ditions different from the standard analysis described in 
Section d] are hsted below. 

1. The ON- and OFF-source data were not matched 
in terms of the geometrical coordinates (azimuth 
and elevation). More data of the longer observa- 
tion times listed in Table [T] without brackets were 
used. The background levels of the ON-source alpha 
distributions were estimated from the OFF-source 
distributions normalized by the observation times 
or from the extrapolation of the ON-source distri- 
butions of alpha > 30°. 

2. In the image cleaning, the threshold for ADC values 
was fixed to be 15 ADC counts. 

3. The box noise cut was not used. Box noise events 
were manually removed looking at image centroid 
and shower rate distributions instead, but no de- 
tails of rejected areas or periods remain. The ADC 
offset correction was not done. The definition of 
the camera edge cut was different and not based 
on the image distance from the camera center. In- 
stead, images for which the centroid was located in 
the area of the outermost pixels were rejected. 

4. The discharge noise cut introduced in Section 15.21 
was not used in the old analysis but instead events 



with centroids located in the bottom quarter of the 
camera were removed since most of the noisy pixels 
were located there. This cut was applied only to 
the 1996 data. 

5. The scaler cut to eliminate the bright star ef- 
fect was more complicated than described in Sec- 
tion 15.21 The mean and standard deviation of 
scaler values with no bright star in the field of the 
pixel were first calculated for each pixel and each 
run. Then signals with scaler values greater than 
8a over the mean value were rejected from the im- 
ages when the distances of the pixel from bright 
stars were smaller than 0?3. 

In the case of the old SN 1006 analysis, calibrated posi- 
tions of the tracking center had been used and therefore 
the accuracy of the object position was better than in 
the case of the old PSR B1706-44 analysis. 

The above analysis procedure has been applied to the 
1996 data, but the box noise cut has been used instead 
of the manual rejection of the noise that is impossible 
to reproduce. The obtained alpha distributions with re- 
spect to the NE rim are shown at the top of Figure [TTl 
where the distribution of OFF-source has been normal- 
ized to that of ON-source using the observation time. 
A broad excess of the ON-source events is found near 
alpha = 0° with the statistical significance of 3.6a be- 
low alpha = 15°. The significance increases to 4.7cr with 
th e different assumption of the background level used 
bv lTanimori et al] (|1998f) . in which the background level 
was estimated from the "flat" region of the alpha plot 
(alpha > 30°) of the ON-source d ata. This result i s sim - 
ilar to the old result obtained bv lTanimori et all ()1998f ) 
in terms of the remaining number of events, the statisti- 
cal significance level, and the alpha peak broadness. The 
middle of Figure [TT] is the alpha distributions of the same 
data obtained with the same method except for the 8a 
scaler cut. The remarkable excess around alpha = 15° in 
this figure has however disappeared further applying the 
discharge noise cut described in Section [5.21 as shown in 
the bottom of Figure [TT] The 8ct scaler cut was origi- 
nally introduced to reduce the effect of the bright stars 
in the field of view, but also has the effect of reducing the 
discharge noise since noisy pixels also have larger scaler 
values. Therefore, it is possible that the excess of the top 
figure consists of discharge noise events remaining after 
the 8ct scaler cut, which is moderately loose for discharge 
noise. 

The above old analysis method has also been applied 
to the 1997 data, and Figure [T^ shows the alpha distribu- 
tions obtained. In contrast to the 1996 result, this 1997 
result is quite different from the old result. No gamma- 
ray-like excess has been found this time and the overall 
features of the old result have not been reproduced. The 
number of remaining background events are about two 
times more than that of the old 1997 alpha plot, and 
therefore some additional selection criteria must have 
been used in the old analysis. However, no further in- 
formation about the old analysis is available. 

6.3. Vela 

It is notable that all of the significances estimated 
in this paper have decreased from those previously re- 
ported. There are two possible reasons for this: the 
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Fig. 11. — Distributions of the orientation parameter alpha for 
the SN 1006 data taken in 1996 obtained using the old analysis 
method. The alpha values have been calculated with respect to 
the NE rim of the SNR. The ON- and OFF-source distributions are 
represented by the histograms with the error bars and the gray solid 
histograms, respectively. Top: distributions obtained using the old 
analysis method described by Kamei (1998), middle: distributions 
obtained using the same method but without the Scr scaler cut, 
bottom: distributions obtained applying the discharge noise cut to 
the middle distributions. 

previous signal for the Vela pulsar region had also been 
enhanced with trials in the parameter domain and/or 
gamma-ray events have been overcut in the present anal- 
ysis since the gamma-ray selection criteria used this time 
are possibly too tight owing to the unmodeled effects in 
the simulations (see the Appendix [A] for more discus- 
sion). The signal from the CANGAROO-I position near 
the Vela pulsar obtained using the 1997 data is enhanced 
from 2.5(7 with the standard analysis in this paper to 
~ 4(7 previously reported by lYoshikoshil (|l998f ) , if the 
size and length selection levels are looser than their stan- 
dard values ( size > 180 ADC c ounts and length < 0?45 
were used bv lYoshikoshil (jl998l )) and more data before 
matching ON- and OFF-source observations (listed with- 
out brackets in Table [T]) are used. We have had no stan- 
dard candle in the southern hemisphere strong enough 
compared to the sensitivity of the 3.8 m telescope and 
cannot fine-tune our simulation code using a real gamma- 
ray signal. It is consequently difficult to distinguish be- 
tween these two reasons. However, in the case of the Vela 
1997 data, 2.5(7 is thought to be a more conservative and 
reliable estimate than the previous value since it has been 
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Fig. 12. — Distributions of the orientation parameter alpha for 
the SN 1006 data taken in 1997 obtained using the old analysis 
method. The alpha values have been calculated with respect to 
the NE rim of the SNR. The ON- and OFF-source distributions are 
represented by the histogram with the error bars and the gray solid 
histogram, respectively. 



obtained with a common analysis method and by reduc- 
ing unaccounted for degrees of freedom. The optimum 
level of the length selection is sensitive to the source lo- 
cation in the camera as seen in Figu re [H and thus the 
length selection level of 0?45 used by lYoshikoshil (|1998l ) 
was conservatively selected considering the search area 
of the 2° X 2° field of view. However, the length selection 
level of 0?45 is very different from the optimized level of 

0?3 for the CANGAROO-I position, and too loose to 
justify the above signal enhancement if we assume our 
simulations to be acceptable for the optimizations of the 
s election levels. 

iDazelev et al.l (|200lD adopted a similar analysis 
method in which they have optimized the parameter se- 
lection criteria so as to maximize the quality factor using 
their simulations. However, they have obtained a null re- 
sult applying their method to the Vela 1997 data. This 
difference is probably due to the different input condi- 
tions in their simul ations which are desc; ribed in the Ap- 
pendix [X] in detail. IDazelev et al.l (|2001h have estimated 
the 3(7 upper limit to the integral flux from the Vela 
pulsar region to be 2.5 x 10~^^ photons cm^^s"^ above 
2.7 TeV, which is about a factor of two smaller than the 
integral flux from HESS J0835~455 obtained by H.E.S.S. 
above the same threshold energy. However, their up- 
per limit is not inconsistent with the H.E.S.S. result 
since only part of the emission region of HESS J0835- 
455 was included in their estimation as in the case of the 
CANGAROO-I result. 

To investigate the above consistency more quantita- 
tively, we have compared the integral fluxes between 
CANGAROO-I and H.E.S.S., changing the integration 
region around the CANGAROO-I position. Figure [13] 
shows the integral fluxes from the Vela pulsar region 
as a function of the radius of the integration region. 
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The circular integration regions have been centered on 
the CANGAROO-I position and Umited to the inside of 
the H.E.S.S. integration region shown in Figure [H] The 
CANGAROO-I fluxes have been calculated considering 
the acceptance distribution in the field of view. The 
source distribution has been assumed to be uniform in 
estimating the average acceptance in the integration re- 
gion, and as shown in Section [5731 the fluxes are reduced 
by about 10% if a Gaussian source distribution is used 
instead. The H.E.S.S. fluxes have simply been obtained 
using the following formula: 

Fr{> iJth) = ^F{> Eth), (9) 

where F{> Eth) and N are the integral flux and excess 
counts for HESS J0835-455, respectively, and Fr{> Eth) 
and Nr are those integrated only inside a radius r from 
the CANGAROO-I position. One can note that the 
CANGAROO-I fluxes exceed the H.E.S.S. fluxes in the 
proximity of the CANGAROO-I position (the most sig- 
nificant excess of 3.0a is found at the radius of 0?2 from 
the 1993-1995 data) but the fluxes are statistically con- 
sistent with each other at large r in the both data sets, 
i.e., the profile of the CANGAROO-I source is more com- 
pact than that of the H.E.S.S. source. This tendency is 
acceptable if the gamma-ray morphology is energy de- 
pendent, since the threshold energies of CANGAROO-I 
are 5-10 times higher than that of H.E.S.S. There is a 
known example, the pulsar wind nebula HESS J1825- 
137, in which the morphology is energy dependent and 
its higher energy e mission is more compact and closer to 
the pulsar position (jAharonian eEaDHoMS)- 

7. CONCLUSIONS 

We have reanalyzed the CANGAROO-I data of 
PSR B1706-44, SN 1006, and Vela using a consistent 
analysis method, in which the gamma-ray selection cri- 
teria have been optimized using gamma-ray simulations 
and OFF-source data. The analysis method is the sim- 
plest using a single cut or two for each image param- 
eter and almost free from arbitrary degrees of freedom. 
The previously reported signals from PSR B1706-44 and 
SN 1006 have not been reproduced with this analysis 
method using the same data, whereas the signal from 
the CANGAROO-I position offset by 0?13 to the south- 
east from the Vela pulsar has still been detected at the 
4.5ct level. The new upper limits to the integral fluxes 
from PSR B1706-44 and SN 1006 at the 95% CL are 
higher than the upper limits obtained by H.E.S.S. and 
so there is formally no inconsistency between them. The 
emission profile of the signal from the Vela pulsar region 
obtained from the CANGAROO-I data has been com- 
pared with that of H.E.S.S. considering the acceptance 
in the field of view of the 3.8 m telescope, which quickly 



drops toward the edge of the imaging camera. The ex- 
pected emission profile for CANGAROO-I obtained by 
reweighting the diffuse emission from HESS J0835-455 
with the CANGAROO-I acceptance distribution in the 
field of view is similar to the observed CANGAROO-I 
sky map, but the latter profile seems to be more compact. 
This difference in shape is possibly due to the fact that 
the threshold energy of the CANGAROO-I observations 
was an order of magnitude higher than that of H.E.S.S. 
and the energy dependent morphology of HESS J0835- 
455 should be studied in more detail. 

We have attempted to reproduce the previous 
CANGAROO-I analyses for PSR B1706-44 and SN 1006 
using the documented analysis methods to investigate 
why the significant gamma-ray signals had been ob- 
tained. "Gamma-ray-like" excesses similar to the old 
results have been found with the analyses carried out 
for the 1993 August data of PSR B1706-44 and the 1996 
data of SN 1006, although they are still less significant 
than the old results. The old results of the other data 
have not been reproduced on the basis of available infor- 
mation, but we can conclude the following. 

1. The previous reports are inconsistent with the re- 
sults obtained here and as some important calibra- 
tion results such as the effect of the ADC offset 
noise were not considered in the former analyses, 
their significance must be questioned. 

2. The previous analyses were unduly complicated, 
resulting in a large number of degrees of free- 
dom in the background event reduction. Changing 
the analysis conditions depending on the data sets 
added to the degrees of freedom, but no reasons for 
the differences were clarified independently of the 
ON-source data. Selecting noise rejection and/or 
gamma-ray selection levels depending on the ON- 
source data itself may have resulted in the appar- 
ently significant signals. 
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APPENDIX 

A. MONTE CARLO SIMULATIONS 

The Monte Carlo simulation code consists of the following three parts: a) generation of EASs, b) Cherenkov light 
emission from the EASs and its at tenuation in the atmosphere, and c) response of the telescope system. To generate 
EASs, GEANT 3.21 (CERNI [19931 ) has been used with the target atmosphere modeled by 80 homogeneous layers, in 
which the densities have been calculated at the average heights using the US standard atmosphere. Primary gamma 
rays and protons have been injected at the height 50 km a.s.l., under which the effect of the homogeneous geomagnetic 
field of 57.7 /iT that has been calculated for Woomera using the International Geomagnetic Reference Field (IGRF) 



Reanalysis of CANGAROO-I data 



17 




Fig. 13. — Integral gamma-ray fluxes from the Vela pulsar region as a function of the radius of the integration region. The circular 
integration regions have been centered on the CANGAROO-I position and limited inside of the H.E.S.S. integration region shown in 
Figure|8] The left and right plots are the results of the 1993-1995 data (Eth = 4.0 TeV) and the 1997 data (_Eth = 2.5 TeV), respectively. 
The CANGAROO-I results are represented by the markers with the error bars, while the H.E.S.S. results are shown by the lines. 

model has been incorporated. Cherenkov photons are approximately emitted from the end point of each charged 
particle path only if they are falling into the telescope mirror area in order to reduce the computation time. Rayleigh 
and Mie scatterings have been considered as the extinction of Cherenkov photons in the atmosphere. The mean 
free path of the Rayleigh scattering is 2974 g cm~^ at 400 nm, and the scale height and m ean free path of the Mie 
scattering have been assumed to be 1.2km and 14km, respectively (jBaltrusaitis et al]ll985f ). Scattered photons have 
been neglected, i.e., deemed to be attenuated. The computation time is further shortened by virtually reducing the 
number of generated photons multiplied by the following factors in advance: 



where N and N' are the numbers of Cherenkov photons before and after the reduction, respectively, and et, e^, e^, and 
Sq are the atmospheric transmission factor due to the above scatterings, the fraction of the azimuthal angle range of 
the telescope viewed by the charged particle, the mirror reflectivity, and the quantum efficiency of the photocathode 
of the PMTs (~ 20%), respectively. 

The telescope optics has been modeled on a parabolic reflector causing the coma aberration on the focal plane 
plus an on-axis Gaussian blur. The blur spot size was measured using a CCD camera, with which an image of an 
on-axis bright star was taken and approximated to the two-dimensional Gaussian of cr = 0?04. Photoelectron signals 
due to Cherenkov radiation and incident on each PMT photocathode are integrated together with the random NSB 
signals and converted to a electronic waveform assuming that the signals have a common triangular pulse shape and 
considering their arrival time differences. Then, event triggers are examined with the conditions of the discrimination 
levels of three photoelectrons for each pixel and five simultaneous hit pixels. The former discrimination level has been 
calibrated comparing scaler values (discrimination rate per 1 ms) exposed to the NSB light with the standard NSB 
brightness 2.55 x 10^'* erg s^^ cm~^ sr~^ (430 ^ 550 nm) compiled by JcUcy (1958). The latter for the pixel multiplicity 
has been determined from lower cutoffs of iVhit distributions as well as its fluctuation of ±1 pixel. One photoelectron 
signal corresponds to 5.0 ±0.7 ADC counts, which has been determined by equating lower cutoffs of ADC distributions 
with the discrimination level of three photoelectrons. This conversion factor is consistent with that calculated from 
parameters of the electronics. Finally, ADC and TDC values of individual pixels are calculated using the 50 ns ADC 
gate width and the 0.24 ns TDC resolution, and recorded in the same format as observed data. 

To check the reliability of the simulations, we have tried to reproduce the trigger rate and image parameter distri- 
butions of background events. The trigger rate R can be calculated using the following formula: 



where / is the differential cosmic ray flux, which is integrated over the simulated energy range between i?min = 500 GeV 
and -Eniax — 100 TeV, Aq = irr^^^cosO is the input area of the maximum radius r„iax = 300 m on the ground with 
the input zenith angle 9, ^Iq = 27r(l — cosG^ax) is the input solid angle of the maximum offset angle Omax — 3° from 
the pointing direction, and e = A^triggcrcd/A^input is the trigger efficiency. The cosmic ray all particle flux has been 



N' = Nete4,ereq, 



(Al) 




(A2) 
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Fig. 14. — Distributions of the image parameters size (top left), width (top middle), length (top right), cone (bottom left), dis (bottom 
middle), and alpha (bottom right) of background hadronic events. The simulated and observed distributions are represented by the 
histograms with the error bars and the solid histograms, respectively. The simulated histograms have been made injecting lO'' protons 
with the power-law energy spectrum of the index —2.7, while the observed histograms have been obtained using the OFF-source data of 
Vela taken in 1997. Simple noise rejections, b < 0.8, size > 200 ADC counts, and Miit > 5, have been applied before filling the histograms, 
which are all normalized to 1 after filling. 



obtained to be f{E) = 3.93 x 10^(£'/1 GeV) ~ ^- ''^ n i~^ sr~^ GeV~^ by fitting a power-law function to the direct 
observation data measured by llchimura et al.1 ()1993D between 500 GeV and 500 TeV. The input primary particles are 
only protons, and so it is possible that the simulated trigger rate is somewhat overestimated in comparison to the 
observed one. The mirror reflectivity and the average zenith angle used here are 75% and 9 — 15°, respectively, to 
compare the result with the Vela 1997 observations. A total of 3775 events resulted in triggers out of 10^ inputs and 
e — 3775/10^ — (3.78±0.06) x 10"'^, where the error is statistical only. Thus, the expected trigger rate is estimated to 
be 4.32 ±0.07 Hz, which is somewhat higher than the real trigger rate of 2-3 Hz. However, this res ult is not inconsistent 
with observations since the fraction of protons in the cosmic ray all particle flux is ab out 40% (Ichimura e t al.l^l993^ 
and the trigger efficiencies for heavier particles are much smaller than that for protons (jAharon ian ct al. 199^. Image 
parameter distributions of the above simulated data are shown in Figure [14] in comparison with those of the Vela 
1997 data. The observed distributions are reasonably reproduc ed by the simulations excep t for the size and cone 
distributions. The same tendency in cone has also been shown bv lDazelev fc PattersonI (|2001h . and they have pointed 
out the possibility that the difference in cone is due to electronic cross-talk. 

Dazclcy & Patterson (2001) also tuned their simulation code by matching trigger rates and image parameter dis- 
tributions between simulations and observations. However, they needed to introduce a signal loss factor to do this, 
which is a free parameter based on the assumption that signal intensities have been lost in the electronics. Using 
such a factor is inconsistent with the above-mentioned calibration results, and our simulations reproduce the observed 
features without it. Low-energy events and the responses to them predominantly determine trigger rates and image 
parameter distributions since the spectrum of primary cosmic rays has a steep power-law form and compensating for 
the disagreement at low energies with a single scaling factor possibly results in distortion at higher energies instead. 
Note that the cosmic ra y spectrum c an be reproduced at relatively high energies using CANGAR OO-I data without the 
loss fa ctor as shown bv I Yoshikoshil ()1999l ). The differences between our simulations and those of iDazeley fc PattersonI 
()2001|) are summarized in Table A major difference other than the signal loss factor is the on-axis blur spot size of 
0?07, which was determined using another CCD measurement. Their larger spot size is possibly due to saturation of 
the CCD, but it is no longer possible to confirm this. Both spot sizes have reproduced shape parameter distributions 
of observed background events reasonably, but they make a significant difference in the narrower width distributions of 
simulated gamma-ray events which are displaced by about 0?02. This is probably the main reason for the differences 
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TABLE 5 

Comparison of simulation details be twee n this work andI Dazeley fc PattersonI 

(I2001I V 



Input Condition This Work Dazelev fc Patterson (2001') 

EAS generation code GEANT 3.21 MOCCA 92 

Composition p 100% p 58%, He 32%, N 10% 

Power-law index of the integral energy spectrum —1.7 —1.65 

Zenith angle at the injection point 15° 25° 

Mirror reflectivity 75% 66% 

On-axis blur spot size'' 0?04 0?07 

NSB rate'' lOMHz/PMT 8.6MHz/PMT 

Signal loss 0% 35% 



^ Standard deviation of the two-dimensional Gaussian. Rate of photoelectron detections with the 
recoated mirror. 

of the optimized gamma-ray selection criteria based on the simulations between this work and iDazelev fc PattersonI 
(|200ll) . 
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